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ABSTRACT 

We examine the proposal that Impulsive x-ray b'.’rsts are produced 
by high-energy electrons streaming from the corona to the chromosphere. 

We find that the currents associated with these streams are so high 
that either the streams do not exist or their current is neutralized 
by a reverse current. 

Analysis of a simple model indicates that the primary electron 
stream leads to the development of an electric field in the ambient 
corona which (a) decelerates the primary beam and (b) produces a 
neutralizing reverse current. It appears that, in some circumstances, 
this electric field could prevent the primary beam from reaching the 
chromosphere. In any case, the electric field acts as an energy exchange 
mechanism, extracting kinetic energy from the primary beam and using it 
to heat the ambient plasma. This heating is typically so rapid that it 
must be expected to have important dynamical consequences. 
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I . Introduction 


Much of the line emission from solar flares must originate in 
plasma of chromospheric density. On the other hand, the corona offers 
a much larger volume for the storage of magnetic energy than does the 
chromosphere. This has led to models in which free magnetic energy, 
released in the corona, propagates to the chromosphere in the form of 
streams of high-energy particles (Sweet, 1969; Sturrock, 1974), Obser- 
vational evidence for high-energy electron streams comes from analysis 
of microwave bursts and impulsive x-ray bursts (Kane, 1972). Evidence 
concerning the origin impulsive x-ray bursts is conflicting and it 
is possible that some radiation originates at the chromosphere, as 
analyzed by Brown (1972), Hudson (1972) and Petrosian (1973 )> and that 
some originates in the corona as proposed by de Jager and Kundu (1963)> 
Kane and Anderson (1970 ) and Kane (I 972 ). 

It has recently occurred to us, and to other solar physicists 
including Brown (1976) and Colgate (I 976 ), that the electrical currents 
associated with the proposed electron streams are so high that the 
streams either do not exist or are neutralized by reverse currents in 
the ambient plasma. The reason for this statement are discussed 
briefly in Section II. An Idealized model in which the primary current 
is neutralized by a reverse current is presented in Section III and 
applied to solar flares in Section IV. Some of the consequences of 
these considerations are discussed in Section V. 
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II. Magnetic Field Associated with Proposed Electron Stream 

If energy released in the corona propagates to the chromosphere in 
the form of high-energy electrons, this electron stream will constitute 
an electrical current. If N is the total number of electrons streaming 
downward over the duration t (sec.) of the impulsive phase, the 
magnitude J (e.m.u.) of the currerit may be estimated from 

J .V. e c“^ N t“^ . (2.1) 

If the transverse and longitudinal dimensions of the stream are 
of order L (cm.), an estimate of the strength B (gauss) of the magnetic 
field produced by the stream is given by 

B 2J l"^ (2.2) 


and the total energy U (erg) of this magnetic field may be estimated 
from 

U ^ B^ ~ L (2.3) 

on 2n 

which becomes 


„ 12-2 .,2 -2 ^ 

U~^ec Nt L 


Kane and Anderson (1970) estimate the total energy involved in a 

29 2 

typical small flare to be 10 ^ srg, the time scale to be ..j 10 sec, 

B 5 

the characteristic length scale to be 10 cm, and infer from the 
x-ray data that the total number of energetic electrons is ~ 10^^. For 


13 2 5 

these values, the above formulas lead to estimates J 10 , B ~ IQ 

34 

and U 10 . It is clear that this model, involving energy transport by 
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an uniieutralized electron beam, leads to unacceptably high values of the 
magnetic field of the beam and magnetic energy associated with the beam, 
When one considers large flares rather than small flares, one faces the 
same difficulty. 

One resolution of this difficulty is to abandon the assumption that 
energy is transported from the corona to the chromosphere by an electron 
stream. One might assume, for instance, that energy is transported as 
heat flux. The above parameters for a small flare then indicate that 
the mean heat flux F (erg cm ^ s would be ~ 10^^. From the formula 
appropriate for a fully ionized plasma (Spitzer, I 962 ), 


F = 


dT 



t.5/2 dT 
dz 


( 2 . 5 ) 


we find that the electron temperature T (°k) must be at least 10^ to 
give the required heat flux. For the concept of heat flux to be valid, 
the mean free path must be small compared with the length scale for 

-1 7.5 

temperature variation, e.g. 10 L, which is ~ 10 for the case now 
being considered. Since the mean free path y for electron-electron 
and electron-ion collisions is given by 


Y 10^ T^ n"^ , (2.6) 

we see that the electron density n (cm must be at least 10^^'^. 

This is substantially higher than is typical in the normal corona at 
the initiation of a flare, although densities of this magnitude do build 
up in the course of a flare. This indicates that, although heat con- 
duction may be an important, and perhaps the dominant, mechanism of 
energy transfer in the later stages of a flare, it cannot be the principal 





mode of energy transfer at flare onset. (One should, however, bear In 
mind the possibility that instabilities may lead to modification of the 
transport properties of the plasma, modifying both the heat conduction 
and effective mean free path.) 


5 


Ill, Model of Electron Stream and Reverse Ctirrent 


It is well known that a plasma tends to preserve charge neutrality. 

A process which tends to give an excess positive or negative charge in 
some region will lead to electric fields acting upon the plasma, Movement 
of electrons in response to this electric field will then restore charge 
neutrality. One expects that analogous processr will also tend to 
maintain current neutrality. If an electron beam is suddenly Introduced 
into a plasma, there is a sudden change in the magnetic field structure 
which will develop induced electric fields opposing the primary current. 
For this reason, we should investigate the possibility that a reverse 
current develops in the corona balancing (or almost balancing) the 
primary current due to an assumed electron stream flowing from the 
corona to the chromosphere. 

Our considerations are substantially simplified if we consider 
a symmetrical flux tube rooted in the photosphere and extending up 
into the corona, and assume that electrons are accelerated at the top 
of the flux tube by the development of stochastic electric fields or 
by some other mechanism. The ejection of these electrons down towards 
the chromosphere then leads to a charge imbalance which leads to the 
development of an electric field of such a magnitude that charge 
balance will be restored. If there is no net change in charge density, 
whatever current flows along the flux tube must be constant along the 
length of the tube. If the tube is symmetrical, this current must be 
zero. Hence the strong tendency of plasma to remain charge neutral 
implies that any primary current will generate a neutralizing secondary 
reverse current. 
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In order to simplify this exploratory analysis, we consider only 
a steady-state situation. We also assume that the background plasma 
can be adequately described by a Maxwellian velocity distribution and 
use transport coefficients based on this assumption (Spitzer, I 96 ?). 

We are interested in the case that the primary electron stream is 
composed of high-energy electrons which consequently have long mean 
free paths. However, we shall find that an electric field develops 
which decelerates the electron stream. But when the electron energy 
becomes comparable with the thermal energy, the mean free path will be 
sufficiently short that the primary electrons will merge with the back- 
ground plasma. As a simple representation of this process, we ignore 
collisions in discussing the primary beam but we assume that an electron 
of the primary beam is absorbed into the background plasma when it is 
decelerated to zero energy. 

If, as a further simplification, we consider a magnetic flux tube 
of constant cross section, we may use the following simple one-dimensional 
form of the Vlasov equation: 


v 5f e ^ 

9s m ds 3v 


(3.1) 


where s measures arc length along the tube, v is velocity (along the tube), 
f(s,v) is the velocity distribution function of the primary electron 
stream, and tp is the electrostatic potential. 

At the top of the loop (s = O), the primary electron stream is 
moving with positive velocity on the side corresponding to increasing 
s. Hence we may, without ambiguity, express f in terms of iji defined by 

P 

mv__ 
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The initial distribution function may therefore be expressed as 


fCO, v) = F(|I() . ■ 

With this initial condition, we find that the solution of the Vlasov 
equation (3-1) 1'^ 

f(s, v) = F(') - cp) . (3.U) 

The current density (e.m.u.) in the primary electron stream 
is given by 

j 

which may be expressed as 

jg = - ^ F(^l( - Cp) d^ . (3.6) 

Since cp will prove to be negative in the region of interest, it 
is convenient to write 





represented by Ohm's law, 




aE 


ds 




where a is the electrical conductivity expreaoed in raodifird gaussian 
units . 

It is convenient to introduce a new independent valrable | to 
replace s, by the relationship 




ds 

cr(s) 


(3.11) 


Then, on substituting equations (3.B) and (3.10) into (3-9) ^nd 
differentiating with respect to we obtain 


^ si = « 

me d ? 


It is convenient to solve this equation for § in terms of 0, 


(3.12) 


§=X(0) , 


(3.13) 


me 

2 


(srs-'"' . 


rather than vice versa. Equation (3.13) becomes 

-2 ^ 
de" 

which may be integrated once to give 

v-l"l /*6 

/ OA 

, “ 1 ria 

'1= 


, 3 . 1 ‘* / 


me 

2 


(SL - g[‘l ■ . 


(3.1'3) 


if we assume that 6 = 0 (tp = O) and X = 0 (? = O) at s .= 0. We find, 
from (3.8), (3.9) and (3.10), that 


Hence equation (3*15) becomes 


IV, Specific Model 


It is now convenient to introduce a specific form for FCijf’t, 




-.1 


This is a power-law distribution at high energy whXch flattens at low 
energy, the "knee" being characterized by 

••O * 1 

We introduce the symbol s) for the flux i>f electrons 'cm s ) 

of energy exceeding eil( at the position s 


s) = I / e/s)) df . 

•'ll) 

If the initial flux is written as we find that 






y + 1 




fio that the total particle flux is given by 


eK ,-v+ 1 


"of'” “ Tfmr *0 

With the form (4.1) for F(\i(), equation (3,17) integrates to ',ive 


X(0) 




(♦o - 


!4.5) 


We easily obtain from (4.5) an expression for the (negative) electric 
potential 6 in terms of the modified distance measure §; 


e(l) = 


^ £_K F 
^0 Y“1 me - 


WY 


- I*- 


0 • 


Hence, from (4.2), we find that 




Y“ 1 i W 


e=K 


^ I '^Q ^ me 

1 1 


Y 


(4.6) 


(4.7) 


i?j,W4NAL PAGE IS P0> 


On noting that the electric current carried by the stream is 
related to H(^, s) by 


jjs) = - f H(0, s) , 




we see that 


Js(s) ’ 


e^K 


{y- 1 ) me 


+ -J- ^ f] ^ 

rO ^ y-l me •M 


(^.9) 


In order to specify the current, particle flux and electric field 

as functions of s, we must adopt a specific form for o(s). A convenient 
» 

approximation to the structure of the solar atmosphere, which is 
expressible in analytic form, is provided by the constant heat flux 
model. If we now assume that s measures distance vertically downwards 
from the corona, and that n = n^, T = T^, at s = 0, this model (Adams 
and Sturrock, 1975) yields the following expressions:' 

1 2/ 7 


T(s) = ~ bFsj' 


n(s) = n^fx^/T(s)) exp j- | | 


■ 10.21 _ 


\ny^ - bFsj^^^ - 


(h . 10 ) 


(h.ll) 


where a « 10 b w lo'^*'^^, and F (erg cm s ) is the downward 

heat flux. 

The electrical conductivity, in modified gaussian units, may be 
derived from the expression given by Spitzer (1962): 


CT = g T 


3/2 


(h.l2) 


-3.6U 

where g « 10 . Hence we find, from (3.11), that | is related to 
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' 4.131 


s by 



Our model is then completely specified by the choice of the coronal 
temperature, the coronal density, the coronal heat flux, Y» the energy 
corresponding to and the injected energetic electron flux. For 
the coronal parameters, we have adopted values typical of the corona 
above an active region (Noyes, 1971): 

T 3 X 10^ K 
n =“ 10^ cm“^ 

F 5 10^ erg cm ^ sec ^ 

We choose to correspond to 25 Icev and we choor,e y = 2.5> corresponding 

to an A' ray spectral index of ~ 3 - 3«5 (Brown, 1975 )> which is a 

reasonable vaiuc; for the spectral index at peak emission (Datlowe, 1975). 

The fraction of uhe beam energy deposited, and the total energy 

6 1 + 

deposited by joule heating between T = 3 X 10 K and T = 3 X 10 K as 

a function of the energetic elecu^'on flux, are displayed in Figure 1, 

For a flare area of 10^^*^ cm^, the uAergetic electron flux inferred 

17 -2 -1 

from a large impulsive x-ray burst corresponds to ~ 10 cm sec 
(Hoyng et al., 1975). Figure 2 illustrates vhe energy deposition 
rate due to joule heating as a function of temperature of the atmosphere 
for this injected energetic electron flux. The or*?inate of Figure 2 
is the time required to raise the ambient plasma temperature by 
10 ^ K. 
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V, Discussion 


For the adopted model, Figure 2 shows that the steady-state assumption 
is not justified for the energetic electron flux required to produce the 
observed x-ray flux by thick target bremsstrahlung. The atmosphere, 
especially the dense low- temperature region, will be heated rapidly to 

T 

a high temperature. Once the plasma is raised to a few times 10* K, 
ohmic losses will become unimportant since the joule heating rate is 
inversely proportional to the electrical conductivity. Joule heating 
is a more effective mechanism for heating the plasma than collislonal 
losses from the energetic electrons, because the ohmic losses are 
caused by thermal electrons in the reverse current which have much 
shorter mean free paths than the energetic electrons. It is likely 
that a larger fraction of the beam energy is deposited than is indicated 

in Figure 1 for two reasons. First, the constant heat flux model under- 

estimates the extent of the low temperature region where the ohmic losses 
are largest. Second, if the drift velocity associated with the reverse 
Current becomes too large, instabilities can develop and reduce the 
electrical conductivity of the plasma. If the conductivity of the 
plasma is low enough, the electric field that develops may be large 

enough to effectively prevent the stream from reaching the chromosphere. 

The time scale for heating the dense low- temperature plasma is 
comparable with the time scales characteristic of impulsive x-ray bursts 
(Hoyng et al., 1975 )- It is possible that thermal bremsstrahlung from 
the rapidly heated plasma can account for a significant portion of the 
observed impulsive x-ray flux. Hence this mechanism offers an explanation 
of the fact that some flares first produce high-energy x-ray emission 



near the top of a loop rather than at the foot-points of the loop 
(Brueckner, I 976 ). Another important consequence of this process is 
that, if thermal emission can account for a substantial fraction of 
the impulsive flux up to ~ 50 kev, then the number of electrons 
required to produce the non- thermal x-ray flux is greatly reduced 
(Brown, 1975). 

The rapid heating of the plasma will also cause a large pressure 
imbalance. The time t (sec.) for the plasma to respond to pressure 
changes by bulk motions can be estimated from 

, (5.1) 

s 

where L is a characteristic length and is the ion thermal velocity. 

'1 2 S 

Even for a temperature as large as 10^ K, this time is long (lO sec) 

compared to the heating time for a length scale of 10^ cm, so that 

the plasma density will not change appreciably during the heating. In 

4 

the constant heat flux model we have adopted, the density at 3 10 K 

6 

is approximately 100 times the density at 3 X 10^ K. If both these 

Y 

regions are heated to a temperature of a few times 10 K, the resulting 
pressure imbalance must be expected to produce explosive evaporation 
of high density plasma into the corona. 

A more realistic time-dependent model is necessary to account 
for the complicated phenomena observed in solar flares. It is also 
possible that collective effects are important and that the assumption 
of a Maxwellian velocity distribution for the ambient plasma must be 
abandoned. However, our analysis of reverse current in terms of an 
idealized model strongly suggests that joule heating is likely to play 

an important role in the heating of solar plasma during a flare. 
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Figure Captions 


.Figure 1, The fraction of the beam energy deposited (solid curve' and 
the total energy deposited (broken curve) by joule heating 
as a function of the energetic particle number flux. 

Figure 2, Time T^j required for ohmic losses to heat the plasma by 

7 

10 K as a function of the ambient plasma temperature. 
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Figure 1 
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FRACTION OF BEAM ENERGY DEPOSITED 
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